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One-particle spectral properties in the normal phase of the two-dimensional 
attractive Hubbard model are investigated in the weak coupling regime 
using the non-selfconsistent T-matrix approximation. The corresponding 
equations are evaluated numerically directly on the real frequency axis. 
For temperatures sufficiently close to the superconducting transition tem- 
perature a pseudogap in the one-particle spectral function is observed, 
which can be assigned to the increasing importance of pair fluctuations. 
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I. INTRODUCTION 

One of the most striking phenomena in high temperature superconductor^! is the 
pseudogap observed in the single-particle spectral density and other quantities way above 
the critical temperature for underdoped samples.il A natural explanation of the pseudogap 
invokes pairing tendencies above the pair condensation scale. The reduced dimensionality 
of the cuprate planes and the short coherence length of Cooper pairs in the cuprate 
superconductors both favor an enhanced influence of pairing fluctuations. 

A popular model allowing the investigation of pairing in two dimensions without com- 
plications coming from other, especially magnetic, fluctuations is the two-dimensional 
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attractive Hubbard mo deli At low temperatures this model has a superconducting phase 
with s-wave pairing. The existence of a pseudogap for single-particle and spin excitations 
in the normal phase of the attractive Hubbard model has been established at intermediate 
to strong coupling by Quantum Monte Carlo (QMC) calculations-BHl For very strong cou- 
pling the pseudogap is a natural consequence of binding of fermions in almost local pairs, 
but the pseudogap has been shown to appear also at more moderate coupling strengths, 
where the momentum distribution function still resembles a Fermi distribution. Pseu- 
dogap behavior in single-particle spectra of the attractive Hubbard model has also been 
obtained within the selfconsistent T-matrix approximation (TMA) at very low densities.! 
In the low density limit two-particle bound states form already at weak coupling in two 
dimensions,0 and the chemical potential can move into the gap between the energy of 
a bound state and the single-particle band. At a band-filling of about ten percent and 
intermediate coupling only some weak tendencies towards gap formation appear within 
the self-consistent TMA,0 and at larger densities standard Fermi liquid behavior seems 
to be restored completely at weak and moderate coupling.00 

In this work we analyze how pairing fluctuations affect the spectral function for single- 
particle excitations in the attractive Hubbard model at weak coupling and intermediate 
(not very small) densities. To this end, we compute the self-energy and the spectral 
function within the non- selfconsistent T-matrix approximation. Using a real-frequency 
formulation of the T-matrix equation^l and exploiting Kramers- Kronig relations we cal- 
culate the resulting spectral function with unprecedented accuracy and energy resolu- 
tion. In this way we are able to detect a dramatic drop of spectral weight at the Fermi 
level for temperatures close to the superconducting transition temperature even at weak 
coupling. This pseudogap behavior is caused exclusively by strong pairing fluctuations 
with small center-of-mass momenta precursing the superconducting instability. It can 
therefore be captured qualitatively within a time- dependent Ginzburg-Landau expansion 
of these fiuctuations.lll The gap is not related to two-particle bound states with large 
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center-of-mass momenta. Our results are in complete agreement with a classical (i.e. 
thermal) pairing fluctuation scenario which has been proposed recently as an explanation 
for pseudogap behavior in the two-dimensional Hubbard model at intermediate coupling 
strength]^ 



II. METHOD 

The Hubbard Hamiltonian is given by 

where and c^^ are the usual creation and annihilation operators for fermions with spin 
a on site and rii^ = cl^c^^^ is the corresponding number operator. The notation 
indicates the restriction of the hopping amplitude t to nearest neighbors on the lattice, 
leading to the familiar form = —2 1 (cos + cos ky) for the free dispersion relation. 
For the attractive Hubbard model the coupling constant U is negative.! 

In this work physical properties are calculated using the non-selfconsistent T-matrix 
approximation (nTMA) described diagrammatically in Fig. |I]. Although standard dia- 
grammatic perturbation theory at flnite temperatures is originally formulated for Mat- 
subara frequencies on the imaginary axis, it is possible to continue the T-matrix equations 
analytically to the real frequency axis before solving them numerically directly for real 
frequencies.i'0 Thereby the difficult and ill-posed problem of analytical continuation of 
numerical results from imaginary to real frequencies is avoided. The retarded real fre- 
quency self-energy within nTMA is then given by 

S(a; + zO+, k) = y y ^ { - /(e) T{e + uj + ^0+, q) ImG°(e + ^0+, q - k) 

+ 6(e) Imr(e + iO+, q) G°(e - - iO+, q - k) } (2) 

with the vertex function 



resulting from the sum over all ladder diagrams with the pair-propagator 

K(. + ,0 .q) = J i^ ,_,..,o_^ + 2^ + ,o^ W 

In the above equations /(e) = (exp(/3e) + denotes the Fermi function, 6(e) = 
(exp(/?e) — 1)^^ the Bose function, and is the propagator of the non-interacting system. 

In principle the Fermi surface is (slightly) deformed by interactions. To avoid the 
computationally expensive self-consistent determination of the interacting Fermi surface, 
we remove Fermi surface shifts by subtracting ReE(0, ki?) from the self-energy E(ct;,k) 
entering the expression for the spectral function, where is a suitable projection of k 
onto the (non-interacting) Fermi surface. The major effect of this is merely a shift of the 
chemical potential. 

The one-particle spectral function is then obtained as 

IN -2S"(cu + ^0+,k) 

A(uj,k) = ^ ^ . (5) 

[u; - (eO - /.) - (S'(c., k) - S'(0, k^))]' + p"(c. + ^0+, k)]^ 

where S' and S" are the real and imaginary parts of the self-energy, respectively. 

We note that eq. (§) is valid only as long as the analytic continuation of the vertex 
function T{z, q) does not have poles away from the real frequency axis. This is discussed 
in more detail in Appendix where equations (^-(^ are derived. It turns out that 
poles in the vertex function appear in the complex plane if the temperature drops below 
a critical value at which the denominator of the vertex function, i.e. 1 — UK{0,0), 
vanishes. This is simply the Thouless-criterion for a superconducting instability.El The 
connection between poles of the vertex function in the complex plane and the onset of 
superconductivity has already been pointed out previously.lli0 For the nTMA the critical 
temperature is identical to the BCS critical temperature One must confirm that 

T > Tc when using eq. (|^) to ensure that the system is in the normal phase. 

Equations (@)-(§) are evaluated numerically. The computation of the imaginary part of 
the pair propagator requires only a one-dimensional numerical integration, since ImG° is 
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proportional to a (5-function. The real part is then obtained from the Kramers- Kronig re- 
lation. The calculation of the imaginary part of the self-energy requires a two-dimensional 
integration, and the real part can again be obtained from the Kramers-Kronig integral. In 
this way the self-energy and spectral function can be computed with very high accuracy 
and resolution. The spectral sum rule is satisfied with an accuracy of the order 10~^ in 
our results. A more detailed description of the numerical treatment of equations (0)-(^) 
is given in Appendix 

For sufficiently large momenta q the vertex function r(ci;, q) has poles on the real fre- 
quency axis, which signal the formation of bound two-particle states. Their contribution 
to the self-energy is investigated separately in order to distinguish between their influ- 
ence and all other effects. In Appendix we show that the domain M"'"^^ in q-space for 
which these bound states exist is given by the condition —2t [cos(ga;/2) + cos(gy/2)] > /i, 
i.e. the region outside a "doubled" Fermi surface. Note that in the case of an isotropic 
dispersion relation the condition for bound states is simply |q| > 2/cir.0 



III. RESULTS 

We now present results for the self-energy and the spectral function. The electron 
density is usually fixed at the intermediate density n = n-^ + rii = 1/2 (quarter-filling), 
unless explicitly stated otherwise. The hopping amplitude is always chosen to be t = 1 
for simplicity, corresponding to a non- interacting band- width Wq = 8. 

We first establish the existence of a pseudogap in the spectral function for a weak 
attraction U = —1.728, for which the corresponding (BCS) transition temperature Tc = 
0.05 is rather low, and discuss its origin. Since we focus on temperatures close to Tc we 
use the reduced temperature = (T — T^/T^. to parametrize T. 

Figure |^ shows the imaginary part of the self-energy k) for a Fermi wave vector 

on the fcj;-axis, k = {kp-, 0) , at various temperatures, with the contribution S^^^. from 
large q two-particle bound states displayed separately for the = 0.01. It can 
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clearly be seen that upon decreasing Tr a sharp peak develops in S" at lo = . This 
leads to a suppression of spectral weight at u = (see below) and thus to deviations 
from Fermi-liquid behavior in the normal phase close to Tc . The origin of this peak is the 
divergence of the vertex function T{uj, q) at = 0, q = in the limit T T^. Stated 
in more physical terms, the peak is thus caused by strong pair fluctuations with small 
center-of-mass momenta in the vicinity of the superconducting instability. 

The bound state contribution ^^bs found to be negligible near u; = and only 
exhibits a rather weak peak at u ^ 2 ^ , which does not affect the spectral properties 
near the Fermi level. This is is contrast to the scenario proposed by Schmitt-Rink et al.,0 
which may be realized only at very low density. The fact that these bound states do not 
generically lead to a breakdown of Fermi liquid behavior has been pointed out already by 
Fresard et al.lli and by Serene.lli 

The resulting spectral function A{uj, k) is shown in Fig. ^ for the wave vector {kp, 0) . 
A pronounced pseudogap forms for sufficiently small reduced temperatures Tr < 0.01. For 
stronger attraction both the width of the pseudogap and the temperature window above Tc 
where the spectral weight at the Fermi level gets suppressed increase, as expected. In Fig. 
^ we show A{uj, k) at k = {kp, 0) for a slightly bigger coupling strength Uc = —2.034, 
corresponding to a critical temperature Tc = 0.1, which is twice as high as before. A 
comparison with Fig. ^ shows that the reduced temperature scale, at which the gap forms, 
is roughly twice as big as in the previous case, i.e. the absolute temperature window for 
pseudogap behavior is four times larger now. The width of the pseudogap increases by 
more than a factor two, upon doubling Tc. The systematic evolution of the pseudogap in 
the spectral function with increasing Tc (i.e. increasing \U\) at fixed = 0.01 is shown 
in Fig. 1^. The width of the pseudogap increases monotonously with increasing \U\, as 
expected, and for fixed it also becomes more pronounced. In Fig. |^ we plot the width of 
the pseudogap, defined as the distance between the two maxima in A{lj, k) , as a function 
of Tc (or \U\), keeping the ratio Tr/Tc fixed at 0.1, such that the pseudogap is equally well 
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developed in all cases. The gap width increases obviously faster than linearly with Tc- 

For other wave vectors on the Fermi surface the behavior does not change much. 
The pseudogap is however less pronounced for k on the diagonal in the Brillouin zone, 
compared to wave vectors on the axes (for fixed parameters). Note that the van Hove 
points are still far from the Fermi surface at quarter-filling. For wave vectors away from the 
Fermi surface the double peak structure in the spectral function disappears and transforms 
into a single peak, as shown in Fig. ^ for the same parameters as in Fig. |^. 

To see how the spectral function changes at lower densities, we now switch to a filling 
factor one tenth {n = 0.2). In Fig. ^ we show results for y4(ti;,k) at k = (A;_f,0) for 
U = —2.667, corresponding to Tc = 0.1 at that density. The qualitative behavior of the 
spectral function is the same as for quarter filling, except for a minute peak with very 
little spectral weight at the lower edge of y4(ti;,k) (see inset of Fig. |]). This peak is 
due to the contribution of poles in the vertex function at large momenta associated with 
two-particle bound states. Satellite peaks in the spectral function caused by bound states 
have also been obtained earlier within the self-consistent T-matrix approximation. i'0 At 
quarter-filling the bound states contribute to the spectral function at an energy which 
lies within the energy range of the much larger contribution from scattering states, and 
is therefore hidden by the latter. 

A pseudogap in the single-particle excitation spectrum of the two-dimensional at- 
tractive Hubbard model at finite (not very small) density has been found previously 
only at intermediate and strong coupling strengths, namely by Quantum Monte Carlo 
calculationsS^ and within a renormalized T-matrix approximation.lli At weak coupling 
the temperature window for pseudogap behavior is very small, and k has to be very close 
to the Fermi surface for a gap to appear in the spectral function A{uj, k). Hence it is very 
hard to resolve the gap at weak coupling in all those numerical techniques for which only 
systems with finite size can be treated in practice. 

Thermal pairing fluctuations have already been invoked by Vilk et al.0 to account 
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for pseudogap behavior at intermediate coupling strengths. These authors showed that 
classical pairing fluctuations can lead to pseudogaps if the correlation length ^ of these 
fluctuations is not only large, but also larger than the thermal de Broglie wave length 
vp/T, where vp is the Fermi velocity. Indeed, we observed at quarter-filling, where the 
Fermi velocity varies already considerably along the Fermi line, that the pseudogap is 
more pronounced at Fermi points with a smaller vp- 

The dominant long-range part of the pairing fluctuations can be described effectively 
using a Ginzburg-Landau expansion of the (diverging) vertex function. Capezzali and 
BeckEi have recently shown that such an expansion captures the singular self-energy con- 
tributions which lead to a pseudogap in the spectral function in two dimensions. 



IV. CONCLUSION 

We have investigated one-particle spectral properties of the two-dimensional attractive 
Hubbard model at weak coupling in the normal phase close to the superconducting transi- 
tion temperature, using the non-selfconsistent T-matrix approximation. For temperatures 
sufficiently close to the transition temperature the increasing influence of pair fluctuations 
induces a pseudogap in the one-particle spectral function for wave vectors near the Fermi 
surface. Two-particle bound states with large momenta do not affect the spectral proper- 
ties near the Fermi level. Our results agree qualitatively with those obtained from Monte 
Carlo simulationJi^ and within a renormalized T-matrix approximationlii for interme- 
diate to strong interaction strengths, and show that a pseudogap also forms for a weak 
attraction in two dimensions. The results show that the general pair fluctuation scenario 
proposed by Vilk et al.i works also at weak coupling. They agree qualitatively with re- 
sults obtained recently within a Ginzburg-Landau treatment of such fluctuations.lli Hence, 
the spectral properties of the two-dimensional attractive Hubbard model seem to evolve 
smoothly when moving from weak to strong coupling, that is there seems to be no critical 
coupling strength at which the behavior changes qualitatively. The size of the energy 



window where pairing fluctuations are important becomes of course very small at weak 
coupling, such that a high numerical accuracy is required to see the effect. In two dimen- 
sions the behavior is thus different than in three-dimensionaS and infinite-dimensionalill 
systems, where gaps open in the normal state only for a sufficiently strong attraction. 

It is clear that the non-selfconsistent T-matrix approximation breaks down in the 
limit T — > Tc, since the propagator is dramatically renormalized at the Fermi level. 
Contributions not contained in the TMA may modify the shape of the pseudogap in the 
limit T Tc, even at weak coupling. The selfconsistent TMA will hardly provide a 
better approximation since vertex corrections are likely to become important as soon as 
self-energy insertions contribute significant ly.lll None of the TMA versions captures the 
Kosterlitz-Thouless physics of the true superconducting phase transition. 

Improved theories should take vertex and self-energy corrections into account on an 
equal footing. An improvement of the TMA which agrees better with QMC results for 
the attractive Hubbard model at intermediate coupling strengths has been proposed very 
recently by Allen and Tremblayii and evaluated in detail by Kyung et al.lii The bare 
interaction is replaced by a renormalized interaction which is determined by consistency 
requirements in this approach. An advantage of this renormalized scheme compared to the 
standard T-matrix is that the temperature window for normal state pseudogap behavior 
becomes larger, in agreement with QMC results at intermediate coupling. However, is 
suppressed too much by fiuctuations in this approach, namely down to zero in two dimen- 
sions, and the expected Kosterlitz-Thouless behavior (at densities n 7^ 1) is not captured. 
The true indeed vanishes only in the half-filled limit {n — > 1), where superconductivity 
would have to break a larger SO (3) symmetry, which leads to stronger order parameter 
fiuctuations and thus to a further enhancement of the pseudogap regime.! By contrast, 
in the standard T-matrix approximation (self-consistent or not) the order parameter fiuc- 
tuations are always underestimated since the phase transition is mean-field like, and the 
pseudogap regime appears thus too small. 
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Functional renormalization group methods, where vertex corrections are automatically 
included, should provide another promising approach to pseudogap physics. Such methods 
have been applied recently to the repulsive Hubbard model,0 but not yet to the attractive 
model. 
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APPENDIX A: DERIVATION OF T-MATRIX EQUATIONS ON THE REAL 

AXIS 

How to continue approximation schemes resulting from a resummation of Feynman 
diagrams to the real frequency axis has been pointed out for example in Ref. |l^, where the 
case of the fluctuation exchange approximation is treated explicitly. The selfconsistent 
T-matrix equations for real frequencies have been presented (without derivation) already 
in Ref. y. Our equations could be obtained from the latter replacing simply G by on 
the right hand side. Nevertheless we will now discuss the analytic continuation of the 
nTMA in some detail, in order to make the paper self-contained, but in particular to 
clarify under which conditions this continuation can be done. 

The non-selfconsistent T-matrix equations for the attractive Hubbard model with 
Matsubara frequencies on the imaginary axis are given by 
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k) = T ^ e^'^-o^ I r(zz/„, q) ^^(zz/^ - q - k) (Al) 

r(^//„,, q) = z ff^. r (A2) 

Kitun, ^) = -Tj2j^ G°(zcu„, k) G°(2Z/„ - lUn, q - k) (A3) 



After carrying out the Matsubara sum in the pair propagator ( |A3D , the analytic contin- 
uation to the real axis (with the correct large frequency behavior) is obtained by simply 
substituting ii/„ i— > z/ + z5, which yields 

^-A ^ M^k /(^-4,,)-/(6°-/z) 
J Atx^ ^ - eq_k + 2/i + 2d) 

The continuation of the vertex function follows immediately as 

The term [/ has been subtracted from F in order to isolate the Hartree term in the 
self-energy 



/J2 
^ (f/ + f(zz/„, q)) G\lVn - ^^n, q - k) 

(A6) 

where the (constant) Hartree contribution can be easily summed to = Un/2. The 
second term in equation (|A6|) is treated using the contour integration technique with the 
integration contour shown in Fig. |^. This yields 

t{iun,\i) = [Tf(0,q)G'°(-zcu„,q-k) 
+ ^ I de { b{e + ii + zw„) f (e + ii + zcj„, q) G°(e + ii , q - k) 
+ h{e + 27') f (e + ii, q) G°(e + ii - iuo^, q - k) 
-(7' ^-7')}! (A7) 
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with 7' > being a small positive number eventually taken to be zero. The term 
Tr(0, q) G'^{—iuJn, q — k) is not covered by any of the contours and is therefore included 
explicitly. 

At this point it is important to notice that for equation (|A7|) to be valid the vertex 
function must not have poles away from the real axis. Such poles appear for \U\ > 
\Uc{T)\ (or equivalently for T < Tc{U)), where Uc is determined by the Thouless criterion 
1 — Uc K{0, 0) = 0. To illustrate this, Fig. [1^ shows the zeros of the imaginary part of the 



pair propagator at z = + i5 for q = together with the real part at these very points. 
The absolute value of the real part has a maximum at = 0. Hence, upon increasing the 
coupling strength a pole first appears for \U\ = \Uc\ at = 0, and then moves into the 
complex plane for larger \U\, along the line defined by the zeros of the imaginary part. 

Provided that \U\ < \Uc\ the steps to finalize the analytic continuation of the self- 
energy can be summarized as follows: (i) use the relation b{z + iun) = —f{z) to replace 
Fermi functions by Bose functions, (ii) make the replacement iujn ^ uo + i'-j with 7 > 7' 
being a small parameter, (iii) take the limit 7' 0''' with 7 remaining finite and use 
limy^o+ 6(e+i7') = —i-n T6{e), (iv) take the limit 7 — * 0+. The result of this procedure 
is 

S(o; + zO+,k) = - [f^ /rfe{ -/(e)f(e + ^ + zO+,q)ImG'°(e + zO+,q-k) 
TT J 47r^ J 

+ 6(e) Imf(e + iO+, q) G\e - to - iO+, q - k) } (A8) 
Inserting T = T — U and adding yields equation (||). 

APPENDIX B: NUMERICAL EVALUATION OF T-MATRIX EQUATIONS 

In order to numerically evaluate the pair propagator the following change of variables 
is made: ^ + ^q-k — —acoskx — b cos ky and A := k^, where a = it cos^q^ / 2) , 
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b = 4tcos(gy/2), and ki = hi — qi/2. Integration in k-space is then written as 



6 Aa 



EQuA 
Qu.l 



di / dX 



d^k . . . = 

a Ai 

where = —(a + 6), ^ = a + 6, and 



6 sin 



(e, A) 



Ai = 0, A2 



Ai = 0, A2 



Ai = arccos 



arccos 



Ao = TT, 



for — a — h<^<a — h 
for a — h <^ <h — a 
for h — a<^<a — h . 



with 1 indicating the summation over all four quadrants of the Brillouin zone. While 
this does not reduce the numerical effort for finite 5 in equation it is a great advantage 
in the limit 5 — 0+ taken in (^). Since the imaginary part of the integrand is proportional 
to 5{ijj + 2/i — the ^-integral is trivial and only a one- dimensional integration needs to 
be performed numerically to compute the imaginary part of the pair propagator: 



A2 



lYaK{uj + zO+,q) 



— TT 



d\ Eg::t(/(/i-4-u)-/(4-^)) 



47r2 



6 sin 



ky A) 



:bi) 



The real part is then obtained using the Kramers- Kronig relation 



ReA> + ,0-. q) = P / ^ lr.K(J + m*,o,) 

J TT Uj' — UJ 

The calculation of the self-energy follows a similar procedure using the slightly different 
transformation^ = — 2 t(cosga; + cosgj,) = e^, A = 2t(cosg2; — cosg^), with the Jacobian 

J(e, A) = 2 (^256 14 + ^4 + ^4 _ 2 (16 t2A2 + le ^2^2 + ' 



after which integration in q-space becomes 

c/2q /■ di 



47^2 



47r2 

U A 



/ d\ J(e, A) ^ 



(B3) 
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with Ai = — (4t — 1^1) and A2 = (4t — |^|). The imaginary part of the self-energy is then 
given by a two-dimensional integral 



ImE(cu + iO+,k) = (B4) 

/ 4^ ^^^^ -li) + h{uj + i-ii)} dX J{^, A) J2 Imr(a; + e - + ^0+, q(^, A) + k) 

and the real part is again obtained using the Kramers-Kronig relation: 

^/ ^ , s n f duj' ImSfct;' + zO^, k) , 
ReE cj + zO+,k = U- + P ^ — — (B5) 

2 . 71 Uj' — U! 



APPENDIX C: TWO-PARTICLE BOUND STATES 

To identify the values q for which two-particle bound states exist we first note that 
the domain, in which the imaginary part of the pair propagator (see Fig. |ll]) is non-zero 
in the limit 6 0+, is bounded from below by uju{q) = — 4t[cos(gx72) + cos(gy/2)] — 2 /i. 
In the case cc;„(q) > the imaginary part jumps at co'„(q) from to a finite negative value. 
This leads to a negative logarithmic singularity in the real part for u c<j„(q)^ and in 
turn to the appearance of a delta-peak in the imaginary part of the vertex function at 
some point z/(q) with < z/(q) < co'u(q), representing a two-particle bound state. The 
border of the region M^^^ in which two-particle bound states exist is defined by ti;„(q) = 
and is shown for quarter-filling in Fig. The condition u;„(q) = reads explicitly 

-2t [cos(g,/2) + cos(gy/2)] = (CI) 

defining a one- dimensional "surface" that is identical to a Fermi surface "up-scaled" by a 
factor two. 
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FIG. 1. Pair propagator, vertex function and self-energy within the T-matrix approximation. 




FIG. 2. Self-energy at U = —1.728 for k = {kF,0) and different temperatures (left), and 
contribution from large q two-particle bound states for = 0.01 (right). 
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FIG. 3. One-particle spectral function A{u;, k) at k = (/cf, 0) for U = -1.728 (Tc = 0.05) and 
different values of the reduced temperature r^. 
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FIG. 4. One-particle spectral function A{u,'k) at k = {kF,0) for U = -2.034 (Tc = 0.1) and 
different values of the reduced temperature r^. 
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FIG. 5. One-particle spectral function A{uj,'k) at k = (kpjO) for different coupling strengths 
(i.e. various Tc) and fixed reduced temperature = 0.01. 
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FIG. 6. Size of the pseudogap as a function of Tc for fixed ratio r^/Tc = 0.1 at k = {kp, 0). 
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FIG. 7. One-particle spectral function A(a;, k) for U 
various wave vectors along the kx-axis. 
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FIG. 8. One-particle spectral function A(a;,k) at k = (A;f,0) for U = -2.667 (Tc = 0.1) and 
different values of the reduced temperature r^; the density is n = 0.2 here. The inset shows the 
tiny peak in the negative frequency tail of A{lo, k) caused by large-q bound states. 
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FIG. 9. Integration contour used to derive equation (|A7|). 
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FIG. 10. Zeros of the imaginary part of the pair propagator at z = to + i6 in the complex 
plane (upper curve) and real part of the pair propagator at these points (lower curve) for q = 
and T = 0.15. 
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FIG. 11. Frequency dependence of the pair propagator for q = (2.6, 2.6) at T = 0.1. 
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FIG. 12. Domain M""^^^ in which two-particle bound states exist (for quarter-fiUing) . 
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